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We study the effect of a potential fourth quark generation on the upper and lower Higgs boson 
mass bounds. This investigation is based on the numerical evaluation of a chirally invariant lattice 
Higgs- Yukawa model emulating the same Higgs-fermion coupling structure as in the Higgs sector 
of the electroweak Standard Model. In particular, the considered model obeys a Ginsparg- Wilson 
version of the underlying SU(2)l x U(l)y symmetry, being a global symmetry here due to the 
neglection of gauge fields in this model. We present our results on the modification of the upper and 
lower Higgs boson mass bounds induced by the presence of a hypothetical very heavy fourth quark 
doublet. Finally, we compare these findings to the standard scenario of three fermion generations. 
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I. INTRODUCTION 



The Sakharov explanation for the matter anti-matter asymmetry of the universe suffers from the 
CP-violating phase of the Standard Model (SM3) falling short by at least 10 orders of magnitude. In 
tIh ■ addition to this concern the Sakharov picture demands a sufficiently strong first order electroweak 

phase transition, which is also objected in the framework of the SM3. However, both of the above 
caveats might be addressable [HHI by the inclusion of a new fourth fermion generation into an extended 
version of the Standard Model (SM4) . Despite the arguments against the existence of a fourth fermion 
generation such a scenario nevertheless remains attractive for two reasons. Firstly, there is a strong 
conceptual interest, since a new fermion generation would need to be very heavy, leading to rather 
large Yukawa coupling constants and thus to potentially strong interactions with the scalar sector of 
the theory. Secondly, it has been argued [l], y[ (and the references therein) that the fourth fermion 
generation is actually not excluded by electroweak precision measurements, thus leaving the potential 
existence of a new fermion generation an open question. 

In our contribution, however, we do not present any statement arguing in favour or disfavour of a 
new fermion generation. Here, we simply assume its existence and focus on the arising consequences 
on the Higgs boson mass spectrum. With the advent of the LHC this question will become of great 
phenomenological interest, since the Higgs boson mass bounds, in particular the lower bound, depend 
significantly on the heaviest fermion mass. Demonstrating this effect will be the main objective of the 
' present work. 

Due to the large Yukawa coupling constants of the fourth fermion generation a non-perturbative 
computation is highly desirable. For this purpose we employ a lattice approach to investigate the 
strong Higgs-fermion interaction. In fact, we follow here the same lattice strategy that has already 
been used in Rcf. for the non-perturbative determination of the upper and lower Higgs boson 

mass bounds in the SM3. This latter approach has the great advantage over the preceding lattice 
studies of Higgs- Yukawa models (see eg. Refs. @-G3 and the references therein) that it is the first 
being based on a consistent formulation of an exact lattice chiral symmetry [14| . which allows to 
emulate the chiral character of the Higgs-fermion coupling structure of the Standard Model on the 
lattice in a conceptually fully controlled manner. The interest in lattice Higgs- Yukawa models has 
therefore recently been renewed • 

The actual details of the considered lattice model will be given in section |TTJ In the following 
section [Hi] the paper elaborates on the pursued simulation strategy for the determination of the upper 
and lower Higgs boson mass bounds, i.e. the selection procedure for the bare lattice parameters. 
The obtained lattice results together with their extrapolations to the infinite volume limit are then 
presented and discussed in section IIVI before the paper ends with some outlook and conclusions in 
section IV] 
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II. THE MODEL 



In order to evaluate the Higgs boson mass bounds we have implemented a lattice model of the 
pure Higgs-fermion sector of the Standard Model. To be more precise, the Lagrangian of the targeted 
Euclidean continuum model we have in mind is given as 

L HY = Pfl + b'0 + -d^d^. + -mfoV + A {ip^f + y b > (?, V) L V>b' R + y v (F, b') L <pt' R 

+ c.c. of Yukawa interactions, (1) 

where we have constrained ourselves to the consideration of the heaviest quark doublet, i.e. the fourth 
generation doublet, which is labeled here (t',b'). This restriction is reasonable, since the dynamics of 
the complex scalar doublet tp (tp = it^p* , 7$ : Pauli-matrices) is dominated by the coupling to the 
heaviest fermions. For the same reason we also neglect any gauge fields in this approach. The quark 
fields nevertheless have a colour index which actually leads to N c = 3 identical copies of the fermion 
doublet appearing in the model. However, for a first exploratory study of the fermionic influence on 
the Higgs boson mass bounds we have set N c to 1 for simplicity. 

Since the Yukawa interaction has a chiral structure, it is important to establish chiral symmetry 
also in the lattice approach. This has been a long-standing obstacle, which was finally found to 
be circumventable by constructing the lattice equivalent of (j) as well as the left- and right-handed 
components of the quark fields t' L R , b' L R on the basis of the Neuberger overlap operator |14l |2~lI [22[ . 
Following the proposition in Ref. [1J| we have constructed a lattice Higgs- Yukawa model with a global 
SU{2) L x U{1) Y symmetry. 

The fields considered in this model are the aforementioned doublet <p as well as N c quark doublets 
represented by eight-component spinors = (f'M 1 6'W), i = l,...,jV c . With D^ ov "> denoting the 
Neuberger overlap operator the fermionic action Sp can be written as 

N c 

s F = J2^ )M ^^ m = v(ov) + dia § (v* > w ) p + + p - dia § (Wf . w ) <t> p - . ( 2 ) 

i=l 

where yt' ,v denote the Yukawa coupling constants and the scalar field <p x has been rewritten here as 
a quaternionic, 2x2 matrix (f>\. — (tp x , fx), with x denoting the site index of the x L t -lattice. The 
left- and right-handed projection operators P± and the modified projectors P± are given as 

p±=^^, p±= ij ^, j 5 = 15 (i-±v^y (3) 

The free Neuberger overlap operator can be written as 

= ptl+-j===y A = VW- P , 0<p<2r (4) 

where p is a free, dimensionless parameter within the specified constraints that determines the radius 
of the circle formed by the eigenvalues of V^ ov ^ in the complex plane. The operator D'W) denotes 
here the Wilson Dirac operator defined as 

where V^ ,&,s are the forward, backward and symmetrized lattice nearest neighbour difference operators 
in direction p, while the so-called Wilson parameter r is chosen here to be r = 1 as usual. 

The overlap operator was proven to be local in a field theoretical sense also in the presence of 
QCD gauge fields at least if the latter fields obey certain smoothness conditions [23], (24|. The locality 
properties were found to depend on the parameter p and the strength of the gauge coupling constant. 
At vanishing gauge coupling the most local operator was shown to be obtained at p = 1. Here, the 
notion 'most local' has to be understood in the sense of the most rapid exponential decrease with 
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the distance \x — y\ of the coupling strength induced by the matrix elements V x °^ between the field 
variables at two remote space-time points x and y. For that reason the setting p = 1 will be adopted 
for the rest of this work. 

The introduced action now obeys an exact global SU(2)l x U(l)y lattice chiral symmetry. For 
f^i € SU(2) and e G IR the action is invariant under the transformation 

t/j ^UyP+^ + Uy^lP-^P, i)^i>p + n\u Y + i)P-U Y , (6) 

<j> ^ Uy<j>Q) L , $^Q, L $U Y (7) 

with the compact notation Uy = exp(zeY) denoting the representations of the global hypercharge 
symmetry group U(l)y for the respective field it is acting on. In the continuum limit Eqs. (J6]{T]) 
eventually recover the (here global) continuum SU(2)^, x U(l)y chiral symmetry. 

Finally, the purely bosonic part S v of the total lattice action S — Sp + S v is given by the usual 
lattice y> 4 -action 

S v = ^Z^lfl^l^x + ^rnl<pl<p x + \(iplip x ) 2 , (8) 

X 

with the bare mass too, the forward difference operator in direction /i, and the bare quartic 
coupling constant A. For the practical lattice implementation, however, a reformulation of Eq. ([5]) 
in terms of the hopping parameter k and the lattice quartic coupling constant A proves to be more 
convenient. It reads 

= -«5>£ [<t> x+fi + + ]T $t $, + A J2 - ^) 2 , (9) 

X , fJj X X 

and is fully equivalent to Eq. ([8]). This alternative formulation opens the possibility of explicitly 
studying the limit A = oo on the lattice. The aforementioned connection can be established through 
a rescaling of the scalar field <& x € 1R 4 and the involved coupling constants according to 

<Pl+i$l\ . A 2 1-2N c X-8k 



W- = V2*( $ S1^ , A=^, <= ± . (10) 



III. SIMULATION STRATEGY 

Due to the triviality of the Higgs sector the targeted Higgs boson mass bounds actually depend 
on the non-removable, intrinsic cutoff parameter A of the considered Higgs- Yukawa theory, which 
can be defined as the inverse lattice spacing, i.e. A — 1/a. To determine these cutoff dependent 
bounds for a given phenomenological scenario, i.e. for given hypothetical masses of the fourth quark 
generation, the strategy is to evaluate the maximal interval of Higgs boson masses attainable within 
the framework of the considered Higgs- Yukawa model being in consistency with this phenomenological 
setup. The free parameters of the model, being the bare scalar mass mo, the bare quartic coupling 
constant A, and the bare Yukawa coupling constants yt>,v thus have to be tuned accordingly. The 
idea for the latter fixation is to employ the assumed fourth generation quark masses mt\h' as well as 
the phenomenological knowledge of the renormalized vacuum expectation value of the scalar field (p 
(vev) . The latter is used to determine the physical scale 1 /a according to 

246 GcV = - = _JL , (11) 



where Zq denotes the Goldstone renormalization constant and v is the lattice vev. 

Concerning the hypothetical masses of the fourth fermion generation quarks, we target here a mass 
degenerate setup with jto# /a = my la — 700 GeV, which is somewhat above its tree-level unitarity 
upper bound [2a ] . The degenerate setting, being unproblematic from a numerical perspective, is owed 
here to the existence of a fluctuating complex phase in the opposite scenario, i.e. in the non-degenerate 
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case [6J. However, it is remarked that we are currently also investigating a set of other mass settings 
to study in particular the quark mass dependence of the Higgs boson mass bounds. 

In lack of an additional matching condition a one-dimensional freedom is left open here, which can be 
parametrized in terms of the quartic coupling constant A. This freedom finally leads to the emergence 
of upper and lower bounds on the Higgs boson mass. As expected from perturbation theory, one also 
finds numerically § that the lightest and heaviest Higgs boson masses are obtained at vanishing and 
infinite bare quartic coupling constant, respectively. The lower mass bound will therefore be obtained 
at A = 0, while A = oo will be adjusted to derive the upper bound. 

However, in the given lattice model the expectation value (ip) would always be identical to zero 
due to the symmetries in Eqs. dSJTJ) . The identification v = (p) would therefore not yield meaningful 
results in Eq. (TTTI) . The reason is that the lattice averages over all ground states of the theory, not only 
over that one which Nature has selected in the broken phase. To study the mechanism of spontaneous 
symmetry breaking nevertheless, one usually introduces an external current J, selecting then only 
one particular ground state. This current is finally removed after taking the thermodynamic limit, 
leading then to the existence of symmetric and broken phases with respect to the order parameter 
(ip) as desired. An alternative approach, which was shown to be equivalent in the thermodynamic 
limit [26l - l28j , is to define the vacuum expectation value v as the expectation value of the rotated field 
ip rot given by a global transformation of the original field ip according to 

pT = U[<p] Vx (12) 

with the SU(2) matrix U[<p] selected for each configuration of field variables {p x } such that 
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(13) 



Here we use this second approach. According to the notation in Eq. ©, which already includes a 
factor 1/2, the relation between the vev v and the expectation value of p rot is then given as 

(0=(J)- (14) 

In this setup the unrenormalized Higgs mode h x and the Goldstone modes g^g^idx^ can then directly 
be read out of the rotated field according to 

(,■"■,:,%.)■ ^ 

The great advantage of this approach is that no limit procedure J — > has to be performed, which 
simplifies the numerical evaluation of the model tremendously. 

The so far lacking prescriptions for calculating the Goldstone renormalization constant Zq, the 
Higgs boson mass m,H, and the quark masses mt>.b> have been discussed in detail in Ref. Q. Here it 
is only stated that Zq is computed from the slope of the inverse Goldstone propagator, the functional 
form of which has been discussed at one-loop order in Ref. [f| . The latter propagator Gq (p) is defined 
as 

Gg{p) = \zZ^v~9%)' (16) 
~9p = —f== e ~ VX 9x (17) 

V L s- L t x 

at the discrete lattice momenta = 2im li /L St t, — 0, . . . ,L s , t — 1. Correspondingly, the Higgs 
boson mass m# is derived from the Higgs propagator given as 

Gh(p) = (h p ~ h - P )i (18) 
h P = . = Ye~ ipx h x . (19) 
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Finally, it is remarked that the lattice results on the quark masses presented in this paper have 
been computed from the decay of the fermionic time correlation functions C / (At) at large Euclidean 
time separations At, where / = t' , b' denotes the quark flavour here. On the lattice the latter time 
correlation functions can be defined as 

Cf(At) = —— g £ £ ( ReTr (f L , t+At , s ■ f R ,t,v) )> (20) 

* s t=0 x,y 

where the left- and right-handed spinors are given through the projection operators according to 

and (F,b')B = (i',b')P- (21) 

It is remarked that the given fermionic correlation function would be identical to zero due to the 
exact lattice chiral symmetry obeyed by the considered Higgs- Yukawa model, if one would not rotate 
the scalar field (p according to Eq. (|12|) as discussed above. This rotation is implicitly assumed in 
the following. Furthermore, it is pointed out that the full all-to-all correlator as defined in Eq. (|20[) 
can be trivially computed. This all-to-all correlator yields very clean signals for the t', b' quark mass 
determination. 



IV. NUMERICAL RESULTS 



For the eventual determination of the cutoff dependent Higgs boson mass bounds several series of 
Monte-Carlo calculations have been performed at different values of A and on different lattice volumes 
as summarized in Tab.fl] In order to tame finite volume effects as well as cutoff effects to an acceptable 
level, we have demanded as a minimal requirement that all particle masses fa = to#, m t > , my in lattice 
units fulfill to < 0.5 and m ■ L s t > 3.5, at least on the largest investigated lattice volumes. Assuming 
the Higgs boson mass m# to be around 500 — 750 GeV this allows to reach cutoff scales between 
1500 GeV and 3500 GeV on a 24 3 x 32-lattice. However, despite this setting strong finite volume 
effects are nevertheless expected induced by the massless Goldstone modes. It is known that these 
finite size effects are proportional to 1 /L 2 S at leading order [2614281 ] . An infinite volume extrapolation 
of the lattice data is therefore mandatory. 



K L a Lt 


N c 


m 2 


A 


Vt' = Vb' 


0.09442 12,16,20,24 32 


1 


2.5910 





3.2122 


0.09485 12,16,20,24 32 


1 


2.5430 





3.2049 


0.09545 12,16,20,24 32 


1 


2.4767 





3.1949 


0.09560 12,16,20,24 32 


1 


2.4603 





3.1923 


0.09605 12,16,20,24 32 


1 


2.4112 





3.1849 


0.21300 12,16,20,24 32 


1 


oo 


oo 


3.3707 


0.21500 12,16,20,24 32 


1 


oo 


oo 


3.3550 


0.22200 12,16,20,24 32 


1 


oo 


oo 


3.1816 


0.22320 12,16,20,24 32 


1 


oo 


oo 


3.1730 


0.22560 12,16,20,24 32 


1 


oo 


00 


3.1561 



TABLE I: The model parameters underlying the lattice calculations performed in this study are presented. 
The setting A = (A = oo) is employed for deriving the lower (upper) Higgs boson mass bound. Depending 
on the lattice volume the available statistics ranges from Nconf = 1, 000 to Nc on f = 20, 000. 

The finite volume data of the renormalized vacuum expectation value v r , the Higgs boson mass 
ran, and the degenerate quark mass ra t > — my resulting from the calculations specified in Tab. U are 
presented in Fig. [TJ These lattice data are plotted versus 1/L 2 and extrapolated to the infinite volume 
limit by means of a linear fit ansatz according to the aforementioned leading order behaviour. Due 
to the observed curvature arising from the non-leading finite volume corrections only those volumes 
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with L s > 16 have been respected by the linear fit procedures. One finds that the intended infinite 
volume extrapolation can indeed reliably be performed thanks to the multitude of investigated lattice 
volumes reaching from 12 3 x 32 to 24 3 x 32-lattices here. 




0.002 0.004 0.006 0.008 0.002 0.004 0.006 0.008 0.002 0.004 0.006 0.008 

U L i U L i 1/^2 



(d) (e) (f) 

FIG. 1: The finite volume data of the renormalized vacuum expectation value v r (a,d), the Higgs boson mass 
run (b,e), and the degenerate quark mass m t / = my (c,f), as obtained from the lattice calculations specified in 
Tab. [I] are plotted versus 1/L 2 a . The upper (lower) row corresponds to the setting A = (A = oo). The infinite 
volume extrapolation is performed by fitting the data to a linear function. Due to the observed curvature 
arising from the non-leading finite volume corrections only those volumes with L s > 16 have been respected 
by the linear fit procedures. 



The quality of the tuning procedure intending to hold the quark masses constant can then be 
examined in Fig. (2b displaying the results of the infinite volume extrapolation of mp . In the considered 
SM4 scenario the fluctuation of the quark mass has been constrained to roughly m t ' = mv = 676 ± 
22 GeV. For later comparisons with the corresponding SM3 scenario we also present the analogous 
summary plot of our earlier investigations 0, HI of the latter setup where the degenerate quark masses 
have been fixed to approximately m t — rrib = 173 ± 3 GeV as demonstrated in Fig. [5^,. 

The infinite volume results of the Higgs boson masses are finally presented in Fig.[3jD. The numerical 
data for the upper mass bound have moreover been fitted with the analytically expected functional 
form of the cutoff dependence of the upper Higgs boson mass bound derived in Ref. [29] . It is given 
as 

up 

^- = A m ■ [log(A 2 /M 2 ) + B m ] - 1 ' 2 , (22) 

with A m , B m denoting the free fit parameters and /i being an arbitrary scale, which we have chosen 
as \i = 1 TeV here. One learns from this presentation that the obtained results are indeed in good 
agreement with the expected logarithmic decline of the upper Higgs boson mass bound with increasing 
cutoff parameter A. 

The reader may want to compare these findings to the upper and lower Higgs boson mass bounds 
previously derived in the SM3. The lattice results corresponding to that setup have been determined 
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FIG. 2: The infinite volume extrapolations of the degenerate quark masses observed in the lattice calculations 
specified in Tab. U are presented versus the cutoff parameter A. In the SM3 scenario (a) the fluctuation of the 
quark mass has been constrained to m t = mi = 173 ± 3 GeV, while m t i = my = 676 ± 22 GeV is adjusted in 
the SM4 scenario (b) . 

in Ref. [U, Q and are summarized in Fig. [5k>. The main finding is that especially the lower bound is 
drastically shifted towards larger values in the presence of the assumed mass-degenerate fourth quark 
doublet. The upper bound is also significantly increased by the fermionic contributions, however 
less strongly. From this analysis it seems conclusive that the usually expected light Higgs boson is 
incompatible with a very heavy fourth fermion generation. 
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FIG. 3: Upper and lower Higgs boson mass bounds are shown for N c — 1, mt = rrib = 173 ± 3 GeV (a) and 
N c = 1, mt/ = my = 676 ± 22 GeV (b). Both upper bounds are each fitted with Eq. (|22|) . The lower bound 
in (a) is also compared to a direct analytical computation depicted by the solid line as discussed in Ref. Q. 
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V. OUTLOOK AND CONCLUSIONS 

The aim of the present work has been the non-perturbative determination of the cutoff dependent 
upper and lower mass bounds of the Higgs boson assuming the potential existence of a very heavy 
fourth quark generation. Due to the potentially strong coupling nature of the associated Yukawa 
coupling constants the lattice approach has been employed here to allow for a non-perturbative in- 
vestigation of a Higgs- Yukawa model. It serves here as a reasonable simplification of the extended 
Standard Model SM4, including a fourth fermion generation. The idea is that the considered model 
contains only those fields and interactions which are most essential for the generation of the Higgs 
boson mass, i.e. the scalar and quark fields and their mutual interactions. The model has been built 
on the basis of Liischer's proposal [l4| for the construction of chirally invariant lattice Higgs- Yukawa 
models adapted, however, to the situation of if being a complex doublet equivalent to one Higgs and 
three Goldstone modes. The resulting chirally invariant lattice Higgs- Yukawa model, based here on 
the Neuberger overlap operator, then obeys a global SU (2)i, x f (1)y symmetry, as desired. 

The fundamental strategy underlying the determination of the cutoff dependent Higgs boson mass 
bounds has then been the numerical evaluation of the maximal interval of Higgs boson masses attain- 
able within the considered Higgs- Yukawa model in consistency with the assumed physical scenario in 
terms of the masses of the fourth quark generation. Owing to the existence of a fluctuating complex 
phase in the non-degenerate case jg|, the t',b' quark masses have been assumed to be degenerate in 
this work. Here we have considered a hypothetical physical scenario with degenerate fourth generation 
quark masses at roughly 700 GeV. To fix the physical scale of the performed lattice calculations we 
have exploited the phenomenologically known value of the renormalized vacuum expectation value 
(vev) of the scalar field 

The aforementioned procedure has actually been simplified in this study by exploiting the - pertur- 
batively expected and numerically confirmed @ - knowledge that the largest attainable Higgs boson 
masses are indeed observed in the case of an infinite bare quartic coupling constant, while the smallest 
masses are obtained in the limit of vanishing bare quartic coupling constant. Consequently, the search 
for the upper and lower Higgs boson mass bounds has been constrained to the bare parameter settings 
A = oo and A = 0, respectively. The resulting finite volume lattice data on the Higgs boson mass, 
the vacuum expectation value of the scalar field, and the quark masses turned out to be sufficiently 
precise to allow for their reliable infinite volume extrapolation. The arising extrapolation results are 
moreover sufficiently precise to resolve their cutoff dependence as demonstrated in Fig. [3] Concerning 
the upper bound the obtained cutoff dependence is in good agreement with the analytically expected 
logarithmic decline and thus with the triviality picture of the Higgs- Yukawa sector. The main find- 
ing, however, is that the upper and lower bounds are both shifted towards larger values. While the 
upper bound is only mildly shifted by roughly 25 %, the lower bound is shifted very drastically to 
roughly 500 GeV. From this analysis it seems conclusive that the usually expected light Higgs boson 
is incompatible with a very heavy fourth fermion generation. 

For the future several connected questions should be further investigated. Firstly, it would be 
interesting to repeat the present study at several other settings of the fourth generation quark masses 
to compute the quark mass dependence of the Higgs boson mass bounds systematically. This should 
be done in particular for the lower bound. Moreover, the implemented model allows to tune the bare 
coupling parameters to arbitrarily large values. It would be of particular interest to determine the 
largest possible quark mass attainable within the Higgs- Yukawa model beyond a perturbative unitarity 
consideration. This can be done by evaluating the model in the limit of infinite bare Yukawa coupling 
constants. 
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